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Abstract 

This expository paper is a tribute to Ekkehart Kroner's results on the 
intrinsic non-Riemannian geometrical nature of a single crystal filled with 
point and/or line defects. A new perspective on this old theory is proposed, 
intended to contribute to the debate around the still open Kroner's question: 
"what are the dynamical variables of our theory?" 

1 Introduction 

In the field of solid state physics, in particular physics of defects, the legacy of Ekkehart 
Kroner who died ten years ago at the age of 81, is invaluable. He has been actively pub- 
lishing for 50 years, mostly as a single author, on the physical understanding of defective 
solids, but also on their mathematical structure. One could make a distinction between 
a first series of paper [17]-[19] where he constructs an original approach to understand 
dislocations, and a later series [20]- [24] where he raises questions, while reporting new 
knowledge in the field. 

Most of the theory can be found in the course [19] but since Kroner also distilled 
many comments, ideas, and computations along other publications, the idea of writing the 
present tribute grew up. It is especially intended to commemorate the 10 th anniversary of 
his death, in order, not to recall (because the author has no privileged relationship with 
Kroner to do so), but to enlighten Kroner's ideas and show how they are found rich enough 
by the author to be diffused, revisited and emphasized today. 

It should be pointed out that Anthony [2, 3] is one of Kroner's direct students who also 
greatly contributed to understand defect lines (in particular, disclinations). Since then, 
many contributions to the field (nonlinear dislocations, dislocation motion, thermodynamic 
of defective crystals etc) have appeared, but surprisingly enough, few only cite Kroner. 
This is probably due to the lack of real school following him, but also due to scientific 
reasons: indeed, Kroner's theory is formulated in physical terms, but appeals to complex 
mathematical concepts, the combination of which is only rarely seen in the literature. It 
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should be emphasized that de Leon, Epstein, Lazar, Maugin and co-authors [11, 14, 8, 25] 
(cf. the well-documented survey [26] and the references therein) have produced significant 
results not only by following, but especially by completing the ideas of Kroner. 

So, the present paper is intended to (i) collect and show Kroner's results in the light of 
a new presentation, (ii) describe the non-Riemannian crystal and show how it can help to 
select appropriate deformation and internal (thermodynamic) variables, (iii) participate to 
the debate around Kroner's question: "what are the dynamical variables of our theory?" 
[22] 

It will be especially stressed that the crystal geometry and the physical laws governing 
defects are inseparable, as is the case in the Einstein's General Theory of Relativity. 
However, we entirely agree with Noll when he writes [27] that "the geometry [must be] 
the natural outcome, not the first assumption, of the theory" (i.e., as in the Continuous 
Distribution of Dislocation (CDD) theory of Bilby at al. [6]). Many geometrical tools and 
mathematical theory required for a rigourous description of the dislocated crystal geometry 
can be found in the landmark papers by Noll [27] and Wang [36], while also pointing out a 
recent book on Continuum Mechanics in that spirit [13]. The approach followed here and 
detailled in [33]- [35] is nonetheless distinct from the CDD theory. 

Single crystals growing from the melt are considered where high temperature gradients 
are unavoidable and hence where point defects are present [32]. Moreover, since there 
are no internal boundaries, the defect lines can take in principle any orientation while 
forming either loops or lines ending at the crystal boundary. However, for the purpose of 
simplicity in the exposition of the theory, we will consider a tridimensional crystal filled 
with a network of rectilinear parallel disclinations and/or dislocations. 

Particular to the chosen approach is the distinction between scales, where the macroscale 
is recovered from the mesoscale by a homogenization process: the singularities (i.e., the 
defect lines) have been erased and hence the density of defects (dislocations and/or discli- 
nations) are recovered by means of smooth fields which we will show responsible for cur- 
vature and torsion of the crystal intrisic geometry. Also, the density of point defects will 
show responsible for the appearance of non-metric terms. In this approach, only objective 
fields are considered to describe defective matter: they are defined across scales although 
their physical meaning might differ. Moreover, no elasto-plastic decomposition and no 
prescription of any reference configuration are required, and there is no assumption of 
static equilibrium (vanishing stress divergence). 

2 Preliminar results at the continuum scale 

Notation 1 In this paper, a scalar, vector or tensor of any order are not typographically 
distinct symbols in the text. The tensor order is specified when equations are written, since 
in this case only, the vector v is written as Vi (with one index), and the tensor U as Uij... 
with a number of indices corresponding to its order. 

The present section focuses on the mesoscopic scale, where dislocations and disclina- 
tions are lines and whose characteristic length is some average distance between neighbor- 
ing defects. The remaining of the medium is a continuum governed by linear elasticity. At 
time t, the body is referred to as 1Z*(t) to represent any random sample corresponding to 
a given crystal growth experiment. In the crystal domain Q,, the meso-scale physics will 
then be represented by a nowhere dense set of defect lines which in 2D are parallel to each 
other. 

Definition 2.1 (2D mesoscopic defect lines) At the meso-scale, a 2D set of disloca- 
tions and/or disclinations C C Q. is a closed set ofQ. (this meaning the intersection with Q. 
of a closed set ofM 3 ) formed by a countable union of parallel lines L^',i elcN, whose 
adherence is itself a countable union of lines and where the linear elastic strain is singular. 
In the sequel, these lines will be assumed as parallel to the z-axis. 
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Since accumulation points (to be understood as clusters of parallel lines) might appear, 
the scale of matter description of this section is named continuum scale. 



2.1 Objective internal fields for the model description 

The present mesoscopic theory is developed from the sole linear elastic strain, which itself 
is defined from the stress field (although the stress-strain relationship is not used in the 
sequel) and therefore is an objective internal field. 

Assumption 1 (2D mesoscopic elastic strain) The linear strain £* is a given sym- 
metric L\ oc (Q,) tensor such that d z £* = 0. Moreover, £* is assumed as compatible on 
Qc := Q\C m the sense that the incompatibility tensor defined by 

INCOMPATIBILITY: n* kl := e kpm ei qn d p d q £^ n , (1) 

where derivation is intended in the distribution sense, vanishes everywhere on Qc- 

In the following definition generalizing the concept of rotation and displacement gra- 
dients to dislocated media, the strain is considered as a distribution on Q (i.e. as acting 
on C\ test-functions with compact support). 

Definition 2.2 (Frank and Burgers tensors) 

FRANK TENSOR: d m ^t ■= e hpq d p £* m (2) 



where xo is a point where displacement and rotation are given. 

Line integration of the Frank and Burgers tensors in Qc (i.e., outside the defect set) 
provides the multivalued rotation and Burgers vector fields u* and 6*. These properties 
are summarized in the following theorem, whose proof is classical. 

Theorem 2.3 (Multiple-valued displacement field) From a symmetric smooth lin- 
ear strain £*j on Qc and a point xo where displacement and rotation are given, a multi- 
valued displacement field u* can be constructed on Qc such that the symmetric part of the 
distortion djU* is the single-valued strain tensor £*j while its skew- symmetric part is the 
multivalued rotation tensor cj*j := —eij k u) k . Moreover, inside Qc the gradient dj of the 
rotation and Burgers fields ui k and b% = u% — eum^*{x m — ii)m) coincides with the Frank 
and Burgers tensors. 

From this result, the Frank and Burgers vectors can be defined as invariants of any 
isolated defect line l/ 1 ' of C 

Definition 2.4 (Frank and Burgers vectors) The Frank vector of the isolated defect 
line is the invariant 



BURGERS TENSOR: 



dibt ■= Ski + £kpq(x P - xop)diu*, 



(3) 




(4) 



while its Burgers vector is the invariant 



B; W := [6J]W = K] W (z) - e klm Qf\x m - x 0m ), 
with [uj k ]^\ [b k ]^ and denoting the jumps of co k ,b k and u k around L 



(5) 
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The three types of 2D defects are the screw and edge dislocation, and the wedge discli- 
nation. As an example, the distributional strain and Frank tensor of an isolated screw 
dislocation (see [33] for the other two) is given by the following L} oc (Q) symmetric tensor 
and Radon measure [1]: 
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Consider a set of countable lines C and remark that the present distributional approach is 
subtle in the sense that the physical condition that be bounded is needed in 

L^ec 

order for d m uj^. to still be a Radon measure [33]. 

Besides their relationship with the multivalued rotation, Burgers and displacement 
fields, the Frank and Burgers tensors can be directly related to the strain incompatibility 
by use of (1), (2) & (3). 



Theorem 2.5 The distributional curls of the Frank and Burgers tensors are 

tiljdldjU3% = V*k 

eujdidjbl = e kpq (x p - x 0p )r]iq, 
with rj* k the incompatibility tensor. 



(6) 
(7) 



From this theorem it results that single- valued rotation and Burgers fields ui* and b* 
can be integrated on Sl if the incompatibility tensor vanishes. 



To complete the model two other objective internal fields are introduced: the disloca- 
tion and disclination densities. 

Definition 2.6 (Defect densities) 

DISCLINATION DENSITY: 8y := ^ tt* (k) T<; k) 6 L(k) (i, j = 1 • • • 3) (8) 

fceicN 

DISLOCATION DENSITY: A*- := V] B* (k) S L(k) (i, j = 1 • • ■ 3) , (9) 

feeicN 

where S L { k ) is used to represent the one- dimensional Hausdorff measure density [1] con- 
centrated on the rectifiable arc with the tangent vector r 4 ' fc ' defined almost everywhere 
on Z/ fe ', while f2*' fc ' and B*^ denote the Frank and Burgers vectors of L^ k \ respectively. 

2.2 Kroner's formula 

In this paper, only a simplified 2D mesoscopic distribution of defects in a tridimensional 
crystal is considered. Accordingly, the vectors 77^, 0£ and A£ denote the tensor components 
Vzki ©zfc an d A* fe . Greek indices are used to denote the values 1, 2 (instead of the Latin 
indices used in 3D to denote the values 1, 2 or 3). Moreover, e a p denotes the permutation 
symbol e za p. 

The contortion as introduced by Kondo [15], Nye [28] and Bilby et al. [6] will show 
a crucial defect density tensor. Kroner [19] understood the importance of this object in 
terms of modelling. 
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Definition 2.7 (2D mesoscopic contortion) 

CONTORTION: K*j := S iz a* - ia£<5y = 1 • • • 3), (10) 

a* := A* — ^at^©*^ — x p)- (11) 

Among several equivalent formulations, the formula relating strain incompatibility to 
the defect densities has been proposed in full generality by Kroner [19], and proved for a 
countable set of 27? lines (this meaning, replacing subscript i by z in the formula) by Van 
Goethem & Dupret [33]: 

KRONER'S FORMULA IN ID: rf k = Q* k + e a pd a Kkp- (12) 

For the expression of incompatibility for a set of skew isolated 3D lines, we refer to [34], 
while general 3D results can be found in [35]. 



3 Preliminar results at the macroscopic scale 

Following Kondo [16], by calling a crystal "perfect", it is meant that the atoms form, in its 
stress-free configuration, a regular pattern proper to the prescribed nature of the matter. 
However, no real crystal is perfect, but rather filled with point and line defects which 
interact mutually. Each defect type is responsible for a particular geometric property, as 
will be described in this paper. In order to reach this crystal, we first need to provide a 
way from passing from the above scale to a scale where the fields have been smoothed. 

3.1 Homogenization 

Homogenization is obtained from the continuum scale by a limit procedure which will not 
be detailled here (cf. [35]), but whose effect is to erase the singularities (isolated ones or 
those resulting from accumulation) and hence to provide a smooth macroscopic crystal. 
Basically we postulate the following limits: 

9*,A*,£* -» 9,A,£, (13) 

where 0, A, £ belong to C°°(fl) and where convergence is intended in the sense of measures 
[1]. The tensor £ will be called macroscopic strain without claiming however that £ is the 
elastic strain (i.e. linearly related to the macroscopic stress a). As a consequence of law 
(13) we directly obtain from (1), (12) and straightforward distribution properties: 

MACROSCOPIC KRONER'S FORMULA: r) k = e a pd a dpUk = Q k + e a pd a Kkp, (14) 

where by (2) and (10), (11), 

MACROSCOPIC FRANK TENSOR: ~8 m Uk ■= e kpq d v £ qm (15) 

MACROSCOPIC CONTORTION: Kij = b~ lz (A* - 6j a eap6 z (xp - £0/3)) - -^A-zSij. (16) 

3.2 External and internal observers 

The external observer analyzes the crystal actual configuration lZ(t) with the Euclidian 
metric gfj 1 = 5ij. The internal observer, in turn, can only count atom steps while moving 
in lZ(t), and parallelly transport a vector along crystallographic lines. According to Kroner 
[20]: "in our universe we are internal observers who do not possess the ability to realize 
external actions on the universe, if there are such actions at all. Here we think of the 



possibility that the universe could be deformed from outside by higher beings. A crystal, 
on the other hand, is an object which certainly can deform from outside. We can also see 
the amount of deformation just by looking inside it, e.g., by means of an electron micro- 
scope. Imagine some crystal being who has just the ability to recognize crystallographic 
directions and to count lattice steps along them. Such an internal observer will not realize 
deformations from outside, and therefore will be in a situation analogous to that of the 
physicist exploring the world. The physicist clearly has the status of an internal observer." 



4 The macroscopic crystal 

At time t, the defective crystal is a tridimensional body denoted by lZ(t). The crystal 
defectiveness is not countable anymore, as was the case in Section 2, since the fields have 
been smoothed by homogenization. However, defectiveness is recovered by the natural 
embedding of the crystal into a specific geometry which will be described in the two 
following sections. 

4.1 Macroscopic strain and contortion as key physical fields 

The macroscopic strain £ and contortion k have been defined by homogenization in Section 
3.1. It turns out that the relevant physical fields are not the Frank and Burgers tensors 
but their completed counterparts [33]: 

Definition 4.1 

COMPLETED FRANK TENSOR QjUJk ■= djU} k — K kj (17) 

COMPLETED BURGERS TENSOR Qjbk : = £kj + e kpq {x p - X 0p )<5jW q . (18) 

The following result is a direct consequence of its mesoscopic counterpart: 



Theorem 4.2 

MACROSCOPIC DISCLINATION DENSITY 
MACROSCOPIC DISLOCATION DENSITY 



Oik = eujdSjUJk (19) 
Aife = € U jdidjb k . (20) 



Vectors n k ,<d k and \ k denote the tensor components n zk ,Q zk and A z k- With the above 
definitions and results, the Frank and Burgers vectors are physical measures of defect 
which are given in terms of the sole strain and contortion tensors: 

Definition 4.3 The Frank and Burgers vectors of surface S are defined as 

MACROSCOPIC FRANK VECTOR. Ofc(S') := / QkdS (21) 

J S 

MACROSCOPIC BURGERS VECTOR B k (S) := / A k dS. (22) 

Js 

As a consequence of (19) and (21) and Stokes theorem, the relation between the completed 
Frank tensor and the rotation gradient appears clear. Moreover, it results from (20) and 
(22) that {5b)j k := djb k appears instead of the displacement gradient. 

In the crystal dislocation-free regions (i.e. where the contortion vanishes), it results 
from the classical integral relation of infinitesimal elasticity that the multiple- valued rota- 
tion and displacement fields read 

u k = u 0k + / S m cj fc d£ m (23) 

J XQ 

rx 

Ui = uoi — EikiLJk(xi — xqi) + / 9 m & fc (£)d£ m . (24) 

J XQ 
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4.2 Bravais metric and nonsymmetric connection as key ge- 
ometrical objects 

A Riemannian metric is a smooth symmetric and positive definite tensor field gij. From 
its symmetry property, there is a smooth transformation a\ such that gij — a™a™<5 m „. 
The metric of the "external observer" on lZ(t) is the Euclidian metric Sij. However, as 
soon as the macroscopic strain £ij is given, another Riemannian metric can be defined on 
lZ(t), namely the 

BRAVAIS METRIC gfj = 8^ - 2£ij, (25) 

where the term "Bravais" (from the notion of Bravais crystal [20]) is to recall that it has 
not a purely elastic meaning. 

The use of this metric on defect-free regions of lZ(t) implies the existence of a one- 
to-one coordinate change between H(t) and IZo, whose deformation gradient writes as 
a-rni = 9mn a ? — Smi — diU m where u m denotes a displacement-like field. Let us remark 
that since small displacements are considered, no distinction is to be made between upper 
and lower indices. 

In the presence of defects, the following object (which is said "of anholonomity" [31]) 
idjk is directly related to the strain incompatibility and hence does not 
vanish as soon as defects are present. This exactly signifies that there is no global system 
of coordinates {xf(xi)} with a smooth transformation matrix a,ji = dixf . In fact, such a 
smooth aij - or, equivalently, such a smooth displacement field only exist in the defect-free 
regions of the crystal. 

Quoting Cartan [9], "the Riemannian space is for us an ensemble of small pieces of 
Euclidian space, lying however to a certain degree amorphously", while Kondo [16] sug- 
gests that "the defective crystal is, by contrast [with respect to the above by him given 
definition of perfect crystal], an aggregation of an immense number of small pieces of 
perfect crystals (i.e. small pieces of the defective crystal brought to their natural state in 
which the atoms are arranged on the regular positions of the perfect crystal) that cannot be 
connected with one other so as to form a finite lump of perfect crystals as an organic unity." 

From the elastic metric, we define the compatible symmetric Christoffel symbols 

BRAVAIS CHRISTOFFEL SYMBOLS rf ;ij = i (digij + djg kl - dugf^j , (26) 

whose torsion r^.ryi := rf ;ij — r^ ;:)i vanishes, while its curvature 

BRAVAIS CURVATURE R?, kmq ~ (d q T^ km + gnp^n;hnJ"f;lq) , (27) 

V / [mq] 

with gn P = S np +£ np the inverse of g„ v under the small strain assumption, and where sym- 
bol [■] denotes the skew symmetric index commutation operator (i.e., A[ mn ] — A mn — A nm ) . 
In the terminology of Wang [36] and Noll [27] a connection such that R^ kmq vanishes is 
called a material connection, while for such a connection ^k-[ij] is denoted as the inhomo- 
geneity tensor. 

Quoting Einstein, "to take into account gravitation, we assume the existence of Rie- 
mannian metrics. But in nature we also have electromagnetic fields, which cannot be 
described by Riemannian metrics. The question arises: How can we add to our Rieman- 
nian spaces in a logically natural way an additional structure that provides all this with a 
uniform character ?" 

In the present case, it is sufficient to replace gravitation by strain, and electromagnetic 
fields by line defects to paraphrase Einstein and raise the question of the apropriate con- 
nection inside the defective crystal. To be complete we should add that in order for the 



7 



theory of dislocations to be closed, it should be combined with the theory of point defects 
which play a role at higher temperature (in the same way as Maxwell theory has to be 
combined with the theory of weak interactions, see Kroner [22]). 

The Bravais geodesies are those lines whose tangent vector n is parallelly transported, 
hence solutions to TjVfri — 0, where V s is the covariant derivative of T B . It turns 
out that on these lines, the internal observer is not be able to recognize any defect line. 
Therefore, the above Bravais connection must be completed by a non-symmetric term. 

The following geometric objects are introduced from the sole dislocation density (or 
equivalently by (18) & (20) from the sole strain and contortion tensors): 



Definition 4.4 

DISLOCATION TORSION: 

CONNECTION CONTORTION: 
NON SYMMETRIC CHRISTOFFEL SYMBOLS: 



Tk-ij '• — Q^ijp^-Jik (^^) 

ATk-ij := Tj.ik + Ti.jk — Tk-ji (29) 
Tk;a := r^--Ar fci y. (30) 



According to Noll [27] Ar fc .[ 3i ] is precisely the crystal inhomogeneity tensor which will be 
shown in the following sections to be directly related to the density of dislocations and 
disclinations. 

5 The macroscopic crystal as a non-Riemannian 
manifold 

By contrast with Kroner's presentation, the present approach shows geometrical objects 
as defined from homogenization of mesoscopic measurable, objective physical fields (20), 
(28) & (29) whose identification with their physical macroscopic counterparts follows as 
(proved) results. 

5.1 Physical and geometrical torsions and contortions 

The following lemma is easy to prove from the definitions. 

Lemma 5.1 The tensor gfj defines a Riemannian metric. The symmetric Christoffel 
symbols Tfij define a symmetric connection compatible with this metric, while T k; ij and 
Ar^ ; ij are skew- symmetric tensors w.r.t. i and j and i and k, respectively. 

The following results makes the link between internal motion of the observer by parallel 
transport and the deformation and defect internal variables encountered. 

Theorem 5.2 (Physical and geometrical torsions) The Cristoffel symbols Tk-,ij de- 
fine a nonsymmetric connection compatible with g t j whose torsion writes as Tk-ij. 

Proof. It is easy to verify [12] that Tk-,ij is a connection since Tf.^j is a connection and 
ATk-ij is a tensor. Denoting by Vfe (resp. Vf ) the covariant gradient w.r.t. T kii j (resp. 
Tk-ij), an d recalling that a connection is compatible with the metric gfj if the covariant 
gradient of gfj w.r.t. this connection vanishes, we find by (30) that 

Vfcflij : = d k gfj -Ti-ikgij -Ti.jkgE = VkgFj + ATi- lk gfj + ATi-jkgu, (31) 

where in the RSH, the 1 st term vanishes by Lemma 5.1 while the 2 nd and 3 rd terms cancel 
each other since ATi-jkgu = ATi-jk = —ATj-^k- It results that the connection torsion, i.e. 
the skew-symmetric part of Ar,-^*, w.r.t. i and k, writes as 

i (AT jlih - AT j]ki ) = -~ {ATi-jk - AT k .,ji) = i((AT*. w - AT t . kj ) + 

(AT k ;ji - AT k;ij ) - (AT i]jk - AT iiW )). (32) 
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Observing that the 1 st term in the RHS side of (32) writes as AT^-ij while, by Definition 
4.4 (Eq. (29)), the LHS and the two remaining terms of the RHS of (32) are equal to 
Tj-ik,Tk-ji and — Tiyfe, respectively, the proof is complete. □ 

Theorem 5.3 (Physical and geometrical contortions) The connection contortion ten- 
sor AFt-ij writes in terms of the dislocation contortion 

ATk-Jj — SkK {6ia6jpe.naKzp) + &ict5 j zCctT K T K + SizSj @E fi T K T k — SkzSiaS j fj£ a fj K Z z ■ 

Proof. For k — z, by Definition 4.4, the last statement of Lemma 5.1, and (16), it is 
found that AF z -ij = ATz-^pSiaSjp, with 

AT z . a f^ — T z . a fj — —-e a fjK z — —-e aT S T fjAz = e aT K T p. 
For k = k, by Definition 4.4 and the last statement of Lemma 5.1, it is found that 
AF K -ij = SicSjp (T K . a /s + Tp- aK + T a -p K ) + 5i a 5j z T z;aK , + 5i z 5jpT z - t p K , 

with T Z -( K = e^ T K, TK and T^ TU = — -e-n/Aj. Since the combination of the terms in Q z vanish 
in AT^-ij, the proof is completed by observing that e Q( 9A K + e Ka A^ = {e a Ke T v)e T i3h. v = 
eaK^fi = e aK K z p. □ 
In conclusion, the non-Riemannian crystal is described from a physical viewpoint by £ 
and k, that is by 15 degrees of freedom. From a geometrical viewpoint the 15 unknowns are 
the 6 components of the symmetric Bravais metric and the 9 (by (28) & (29)) nonvanishing 
components of the connection contortion. 



5.2 The Bravais crystal 

The following definition introduces two differential forms whose path integrations gener- 
alize (23) & (24) to the defective regions of the crystal. 

Definition 5.4 (Bravais forms) 

dujj := "dpujjdxp, (33) 
dPu := -Tt.^dxf}. (34) 

In the literature the existence of an elastic macroscopic distortion field is generally postu- 
lated together with the global distortion decomposition in elastic and plastic parts (for a 
rigorous justification of the latter, see [23]). The present approach renders however pos- 
sible to avoid this a-priori decomposition. Nevertheless, the following theorem introduces 
rotation and distortion fields (which must not be identified with the rotation and distortion 
as related to the macroscopic strain) in the absence of disclinations. As a consequence and 
in contrast with the classical literature where it is basically postulated that dislocation 
density is the distortion curl, this relationship is here proved. 

Theorem 5.5 If the macroscopic disclmation density vanishes, there exists rotation and 
distortion fields defined as 

BRAVAIS ROTATION OJj[x) := (J° + / dojj, (35) 

J XQ 

f' X 

BRAVAIS DISTORTION f3 ki (x) := £ k i(x°) - e kh uj" + / d/3 kh (36) 

J x 

= £ k i(x)-e k ijUj(x) (37) 

where ujj is arbitrary and the integration is made on any line with endpoints xo and x. 
Moreover, 

dafikp — d a £k0 + e kp pd a u) p and e a pd a fikf) — A. fc . (38) 
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Proof. By Definition 4.4, the symmetric part of the connection writes as 

- ^(i-,k)pdxp = - ^dpg^dxp = -^d m g k idx m = d m £ k idx m = d£ kl , 
while, by Definition 4.4 and Theorem 5.3, the skew-symmetric part writes as 

Observing, by (17) and Definition 5.4 and Theorem 5.3, that dojj = Qpujjdxp = —^tikjF^.^pdxp, 
it results that d/3ki = d£ki — tkljdwj. Under the assumption of a vanishing macroscopic 
disclination density, the existence of single-valued Bravais rotation and distortion fields 
follows from (33), (19) & (21). Moreover, since dafikp = d a £k@ — ekpjdaUj, Eq. (38) is 
satisfied by (19) & (20). □ 

Remark 1 Eq. (35) indicates that symbol in (33) becomes a true derivation operator 
in the absence of disclinations. 

Remark 2 Referring to "Bravais" instead of "elastic" rotation and distortion fields is 
devoted to highlight that these quantities do not have a purely elastic meaning. In fact, the 
Bravais metric is not even needed since the internal observer only requires the prescription 
of the connection, and subsequent path integration of the forms: 

d£ki = — T^.^pdxfj, dojj :— — -tikj^ [i-,k]pdx p and dfikl := —Ti^pdxp. 

Remark 3 The Bravais distortion does not derive from any "Bravais displacement" in 
the presence of dislocations. In fact, around a closed loop C , even if the disclination density 
vanishes, the differential of the displacementas du k := /3 ka dx a verifies by Theorem 5.5: 

/ du k = I e a pd a PktidS = / Qpbkdxp = A fe (S) = / e a pd a Qpb k dS. (39) 
Jc Js Jc Js 

Remark 4 Theorem 5. 2 defines an operation of parallel displacement according to the 
Bravais lattice geometry. The parallel displacement of any vector Vi along a curve of 
tangent vector dx& is such that dx^VaVi = and hence that the components of Vi vary 
according to the law d^Vi = —Ti-jpVjdx^ [12]. This shows the macroscopic Burgers 
vector and dislocation density together with the Bravais rotation and distortion fields as 
reminiscences of the defective crystal properties at the atomic, mesoscopic and continuum 
scales. In fact, ifdxv,dx^ are two infinitesimal vectors with the associated area dS := 

tv^dx^ dx^\ it results from Eq. (28) and the skew symmetry ofF k;a p that, in the absence 
of disclinations, 

dB k = A za dS = -e a pT k .i} a dS = -T k; p a (dx^dxy - dx^dx^P), 
whose right-hand side appears as a commutator verifying the relation 
dB k = Capda^pbkdS = -e a ^d (a) (dx (,3) ). 

5.3 Motion of the internal observer 

The internal observer will be represented by the k th geodesic basis element e^{x) solution 
to 

(ek)jVj(ek); = (with no summation on k), (40) 

where V is the covariant derivative of Y as given by (30). We have seen in the above 
two sections that it was sufficient to provide him with a connection, i.e. with a law of 
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parallel transport inside the crystal. In fact at this stage the internal observer is not able 
to measure distances, while he can measure the disclination (resp. dislocation) content of 
a surface S by boundary measurements of O^fefc (resp. Qpuik) on the curve C enclosing S 
(which depend merely on F - cf Remarks 2-4). 

The notion of metric connection can be explained as follows. Let the external observer 
be equipped with the Bravais metric and Cartesian coordinate system {a;;}. Since Ti-,k m = 
V m (ek)a, we have on a portion A — B of geodesic k, 

(e k )i(B) - (e k )i(A) = / Ti- km ,dx m = lim r !;tm (i')(e k ) m (a;')As ! , 

A l<i<N 

where x 1 are discretisation points on the curve with endpoints x 1 = A and x N = B, and 
As 1 a tending to zero element of the geodesic. Moreover, if the connection is compatible 
with the metric g B , the angles between these lattice vectors and their (unit) length remain 
invariant during parallel transport. So, we understand Kroner [21] when he says: "when a 
lattice vector is parallelly displaced using F along itself, say 1000 times, then its start [say, 
A] and goal [say, B] are separated by 1000 atomic spacings, as measured by g B . Because 
the result of the measurement by parallel displacement and by counting lattice steps is the 
same, we say that the space is metric with respect to the connection r." 

Moreover, as long as (e k )i(A) (the internal observer) is parallely transported along a 
closed curve C with start- and endpoint A, the gap as created when he comes back to his 
orgin can be measured by the external observer, since by Stokes theorem 

( e k")z — ( e k)i = / ^l-.kmdXm = / e pqm d q T ' 1-kmdSp 

Jc Js 

tpqm (VqTl ; km + (Tl;p m Tp;kq + r 

I \pq± p\km 

+ F 

) dS p , 

where (ek")a denotes the base (ei)fe after being parallely transported along C. Since the 
term inside the parenthesis is symmetric in m and q, we have [12]: 

= / tpgm — (V[ 9 V m j (ek)i + \7 p(ek)iTp; m q) dS p = I ep qm — Ri- nmq (ei t .)ndSp, (41) 
Js 2 Js 1 

with the definition of the Riemannian curvature tensor 

Rl;kmq '■= R^kmq + ^Rl;kmq, (42) 

where by (27) & (30), R B and AR denote the Riemann curvature tensors associated to 
r s and AF, respectively. 

By (41), the internal observer is convinced to return to his startpoint while the external 
observer however can see the gap as created by the crystal curvature, itself resulting from 
the presence of defects. 

5.4 Geometric and physical curvatures 

Let us remark that in the absence dislocations (T — AR — A = 0), the gap is merely 
due to curvature effects with a curvature tensor directly related by (14) to the disclination 
density by Ri- kmq = —eikiemqjOij ■ It should however be noted that in the absence of 
disclinations, the curvature is not vanishing but depends on the sole contortion, since from 
(14) & (42), 

kmq — tlki^mqj ^ipn 

dpKj n + ARl-kmq, (43) 
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where by Theorem 5.3, AR is linearly related to the contortion. It is computed from (43) 
that the Ricci and Gauss curvatures [12] read 



RICCI CURVATURE 



GAUSS CURVATURE 




(44) 



(45) 



while Einstein tensor reads 




? — R., 




(46) 



in the presence of dislocations and disclinations, thereby contradicting Kroner who iden- 
tified Einstein tensor with the disclination density in [21]. 

Moreover, since the macroscopic strain can be decomposed into (symmetric) compat- 
ible and (symmetric) solenoidal parts [33], where only the second one as denoted by £ a 
is relevant for the incompatibility tensor, it results that its trace £p P satisfies by (45) 
— AS" — R B , thereby showing how the Gauss curvature is related to the variation of 
matter density. 

5.5 Summary of the non-Riemannian metric crystal 

The crystal equiped with {g B , T} has the following properties: (i) the geodesies of F are 
the crystallographic lines; (ii) the effect of parallel displacement of the internal observer 
(equiped with F) along a crystallographic line is equivalent to counting the lattice steps; 
(iii) the defect content, i.e. disclination and/or dislocation densities can be computed 
from measures of T only; (iv) the torsion of V is merely due to the presence of dislocations, 
while its curvature is due to the presence of both disclinations and dislocations; (v) in the 
absence of disclinations, there exists a single-valued rotation and distortion field; (vi) if 
and only if there are no defect lines, T is Euclidean and there exists a holonomic coordi- 
nate system. In the latter case only, one can properly speak of a reference configuration, 
of single-valued rotation, displacement and distortion fields, with the macroscopic strain 
compatible with the displacement field. 

Figure 1 illustrate the inseparable link between physics and geometry. On the one 
hand, the physical fields can be set apart: the deformation and defect internal variables 
are shown in rectangular and hexagonal boxes, respectively. On the other hand, the purely 
geometrical object are in oval boxes. The double boundary line means that the quantity 
contains differential combinations of other fields (as connected by arrows), while single 
lines mean algebraic combinations only. 

The main deformation field is the strain, while the distortion and rotation are only 
obtained in the absence of disclinations (see Theorem 5.5). Because they depend on an 
arbitrary point where there value is assumed known, they are considered inappropriate 
as model variables. Concerning defect internal variables, one could indifferently chose the 
dislocation torsion or contortion. Let us mention that since strain instead of distortion is 
chosen, the deformation and defect variables should be considered as independent physical 



5.6 Nonmetricity, teleparallelism and the paradox of the fiat 
crystal 

Let us first remark that the notions of metric and of connection must be considered as 
distinct. This has been emphasized in [7] where it is recalled that historically it has not 
been so for a long time (including Einstein literature). Here, we have seen that the metric 
is attached to the notion of external observer, while the connection is attached to the 



fields. 
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Figure 1: Link between physics and geometry of defects 



notion of parallel displacement of the internal observer inside the crystal. We have seen 
that parallel displacement with F B and counting-step measurements are the same in a 
crystal filled with line defects. This was true because F B was compatible with g B and 
hence crystallographic basis elements remain crystallographic as transported along the 
crystallographic lines. 

Suppose now that the crystal also contains point defects. According to Kroner [23], 
"nonmetricity means that length measurements are disturbed. It is easy to see that this 
just occurs in the presence of point defects. In fact, when counting atomic steps along 
crystallographic lines to measure distance between two atoms, [the internal observer] feels 
disturbed when suddently a vacancy or an intersticial emerges instead of another atoms 
[of the perfect crystal] ." 

Let Cv and Ci be the scalar vacancy (resp. interstitial) concentration, that is the 
number of vacancies (resp. interstitials) per unit volume of crystal. Then, the following 
metric as proposed by Kroner [20]: 

g' = (l + Ci-Cv) 2 g B , (47) 

verifies dV = \/det g'dVo = (1 + AC)\/det g B dVo, with dVo the volume element of the 
stress-free crystal and dV that of the actual one, and where AC = Ci — Cv is the excess 
atomic content of dV . An evolution equation for Cv and Ci (and hence for AC and g') 
will be given in Section 6. 

It is clear that the non-metricity defined as Qj-i k := Vjg'ik 7^ [31, 21] must now enter 
the geometric point- and line-defect model. Differentiation V is here intended with respect 
to connection F, as defined by 

CHRISTOFFEL SYMBOLS WITH POINT DEFECTS: F k ;ij := F' k . i:j — AF k ;ij — ^SF k ;ij (48) 

with Ar fe; ij given by (28) & (29) and where F' k . i:j := \ (dig' kj + djg' ki - d k g'ij) and 
NONMETRIC CONTORTION: 5F k -ij := Q r<lk + Q iijk - Qk;ji- (49) 
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Since Q is a tensor quantity, it is expected to play a role for the physical description of 
the crystal (and to obey an evolution equation as related to the other defects). 

The paradox of the flat crystal is the fact that a defective solids with in addition to 
defect lines a certain amount of point defects can recover a vanishing curvature if the 
following balance holds: 

TELEPARALLELISM ASSUMPTION: SRi-kmq = -(R'l;kmq + ARi ;kmq ), (50) 

where the three terms are the curvature of 8F, F' and Ar, respectively. This is what 
Kroner (and Bilby et al [6]) calls teleparallelism, by this meaning that the global connection 
curvature vanishes and hence that the internal observer ends up parallel when travelling 
along a loop. It should be emphasized that teleparallelism is often considered as a working 
assumption [6, 23]. 

However, we rather follow Kroner [20] when he says that "curved crystals are possible 
only if the curvature is, in some sense, compatible with the considered crystal structure" , 
which means that instead of a flat crystal given by (50), the connection curvature R of the 
actual crystal should be such that point and line defects accomodate to satisfy 



Particularizing (51) we learn from identity R(i. k ) mq + V [m Q q] . lk + T p . mq Q p -i k = [31, 21] 
that point defects and dislocations must be geometrically related, which phenomena is 
well known from solid-state physicist: "dislocations moving perpendicular to their Burgers 
vector produce point defects, and similar processes occur when dislocations cut each other" 
[20] (concerning non metricity, see also [4]). 

6 Concluding remarks: the choice of model vari- 
ables 

Let us conclude by attempting to answer Kroner's question of the Introduction. To the 
knowledge of the author, this question has not been answered yet, or to say the least, there 
is still no agreement on the answer. 

In fact it depends of the physics which one wants to capture. If motion of dislocations 
is modelled, then not only the conservative glide but also the non-conservative climb 
mode must be taken into account. Non conservation is due to the presence, creation, 
anihilation and motion of point defects, and these processes require high temperatures and 
non-negligible temperature gradients [32, 30]. Therefore a complete model of dislocation 
motion must be thermodynamical, away from thermal equilibrium (i.e. irreversible), and 
coupled with the motion of point defects. In [32] the following set of PDEs showed very 
good results to model point defects: 



with the Lagrangian derivative D/Dt, and where Ck, Dk, Dk and P mean (scalar) con- 
centration, (tensorial) equilibrium diffusion, thermodiffusion and (scalar) recombination 
(K — I for interstitials and K = V for vacancies). 

Concerning the motion of dislocations (we assume here that disclinations are negligi- 
ble), a similar equation as (52) should be proposed, with an inter-dislocation recombination 
term and a term of interaction with point defects, collectively denoted by P (also appearing 
in (52)). However, the dislocation density cannot be scalar (which is the case in most of 
the current models available in the literature) but must be the tensorial A (or equivalently 
the contortion k). The PDE could read 




(51) 



DC K 
Dt 



V • [DkVCk + DkCkVT^ - P, 



(52) 




(53) 
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with appropriate boundary condition and where D and D are tensor diffusivities of order 
4. Moreover the contortion verifies the conservation law [33] V • k = V (tr k), meaning 
that the mesoscopic dislocations are loops or end at the crystal boundary, in such a way 
that Eq. (53) amounts to a system of 6 coupled PDEs. Moreover the expression of P must 
somehow satisfy the geometric interaction between defects as given by (51). 

Concerning the deformation variables, let us first observe Figure 2. Incompatibility is 



the final quantity as obtained by recursive differentiation of either the strain (twice), or 
the contortion (once). It hence shows the ultimate convergence of the initially set apart 
deformation and defect variables. The other oval boxes denote defect variables obtained 
from the two key model variables: strain and contortion. 

It should be observed that all relations between strain £, Frank tensor dui, contortion 
k and incompatibility are obtained by means of recursive application of the curl operator 
(either Vx or xV). 

As a first step, Kroner proposed an ("athermal") Gibbs free energy reading W = 
W(F e ,A) with F e the elastic deformation gradient. Restricting hence to statics, he 
thereby attempted to answer his question [22]: "what are the independent (extensive) 
[explicit state] variables entering the free energy (at constant temperature)?" However in 
[23], he recognizes that the use of F e is inevitably ambiguous because the elasto-plastic 
decomposition is not unique. 

According to our theory, the free energy naturally reads from the diagram on Fig. 2 as 
a first strain gradient (see [10, 25]) model W — Wi(£,doj; n) where the strain gradient is 
however replaced by its curl. Equivalently it could read W = Wi{£ , Qoj, A) by combination 
of the last two variables of Wi, or even W = Wz{£ , 86, A) by combination of the first two 
variables of W%. Let us observe that according to Theorem 5.5, W3 can nevertheless 
be compared with W as soon as 36 is identified with a distortion (i.e. a deformation 
gradient), although not becessarily the elastic one (cf. Remark 2). Moreover, a curl 
differential relation is also observed between W3 last two variables (cf. Remark 3). A 
recent thermodynamic analysis with W3 has been remarkably reported by Berdichevsky 
[5] where 9b is identified with the plastic distortion. Let us remark however that by (17) 
& (18), 86 is not, because of the prescription of the arbirary xq, an unambiguous state 




Figure 2: The deformation and defect state variables 



15 



variable, as opposed to duj. 

They are however reasons to be tempted by the choice W = Wi because (i) all variables 
are explicit state variables defined by objective fields (which do not appeal to reference 
configurations, arbitrary plastic parts, or points xo), (ii) there is a distinction between 
(strain-like) deformation and (internal) defect variables, (iii) all variables have clear and 
unambiguous physical and a geometrical meanings. If needed, all other variables (such as 
V£ , Set), 36, A, 77) can be recovered as implicit state variables of the model [22]. Moreover, 
if applying the curl operator twice to the strain and once to the contortion, then the only 
additional model variable naturally appearing is the incompatibility, of both deformation 
and defect nature. This sounds like a closure on the recursive iteration for (higher-order) 
models. 

Nonetheless, we rather prefer to introduce incompatibility through Kroner's formula 
(14) as a constraint to the Gibbs energy with 6 degrees of freedom as coupling strain, 
Frank tensor and contortion: 

Vx^xV = Vx9k) = cxV, (54) 

while the equilibrium law writes as the following equation with 3 d.o.f.: 

dW 

V with cr:=^" and W = W(T, VT; £, dui; k), (55) 

where / is the sum of external forces and of configurational (internal) pseudo-forces di- 
rectly related to k and to the derivative of the so-called dislocation moment stress dWi /dn 
[23, 5]. Let us also mention that the additional constraint of incompressibility must be 
added in order to avoid climb and point defects [22]. Also, a remarkable discussion on the 
nature of W = Wi(db) (with identification of 9b with the distortion) in a nonlinear and 
variational setting can be found in [29]. 

Summarizing, an athermal model of dislocations requires to solve equations (53)- (55) 
which involve a total of 15 degrees of feedom. This number is exactly the number of 

d. o.f. required by the internal observer to parallel displace inside the crystal (through 
the nonsymmetric connexion). Moreover to be closed the theory must involve point and 
line defects, and hence must consider high temperature and temperature gradients. So, as 
recognized by Kroner [23], the dislocation model must be cast within the general frame 
of irreversible thermodynamics because the time variations of the internal variables create 
thermal dissipation. 

This is the main reason why a huge work remains to be done in order to determine, 

e. g., the stress-strain relation, all other constitutive laws, the diffusion coefficients (which 
depend on the crystal internal symmetries, glide planes, etc), and the defect interac- 
tion/production terms. 

The author is unable to answer definitely any of these questions but will pursue research 
on the topic. This paper aims at recalling Ekkehart Kroner's legacy, and in particular the 
fundamental questions he raised which are still open and crucial nowadays. It is also 
aimed at stressing that solutions to dislocation modelling will most probably arise from a 
strong interplay between mathematics and physics, as remarkably done by Kroner along 
his papers [17]-[24]. 
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